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ABSTRACT 

We present a new set of dissipationless N-body simulations to examine the feasibility of 
creating bright ellipticals (following the Kormendy relation) by hierarchically merging 
present-day early- type dwarf galaxies, and to study how the encounter parameters 
affect the location of the end-product in the (/Xe) — Re plane. We investigate the 
merging of one-component galaxies of both equal and different masses, the merging of 
two-component galaxy models to explore the effect of dark halos on the final galaxy 
characteristics, and the merging of ultracompact dwarf galaxies. We find that the 
increase of (fie) with i?e is attributable to an increase in the initial orbital energy. 
The merger remnants shift down in the (/ig) — Re plane and fail to reach the KR. 
Thus, the KR is not reproducable by mergers of dwarf early-type systems, rendering 
untenable the theory that present-day dwarfs are responsible for even a small fraction 
of the present-day ellipticals, unless a considerable amount of dissipation is invoked. 
However, we do find that present-day dwarfs can be formed by the merger of ultra- 
compact dwarfs. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: dwarf - galaxies: for- 
mation - galaxies: fundamantal parameters - methods: numerical 



1 INTRODUCTION 

Despite the progress made in understanding the physics 
of early- type galaxies, a single widely accepted formation 
model, capable of explaining all their observed properties, 
does not yet exist. There is a variety of scenarios describing 
how these objects were formed, ranging from the collapse of 
clumpy protogalaxies (Eggen, Lynden-Bell & Sandage 1962, 
Larson 1975) to the merging of smaller galaxies (White & 
Rees 1978, Kauffmann 1996, Cole et al. 2000). Nevertheless, 
the remarkable fact is that early-type galaxies demonstrate 
a very tight kinematic-structural relationship, usually re- 
ferred to as the Fundamental Plane (Djorgovski & Davis 
1987, Dressier et al. 1987), and any theory of galaxy forma- 
tion and evolution must be able to account for its tightness 
and environmental independence (de la Rosa, de Carvalho, 
& Zepf 2001). 

The Fundamental Plane (FP) combines photometric pa- 
rameters {Re, the effective radius, along with {ne), the mean 
surface brightness within R^) with a spectroscopic observ- 
able (the line-of-sight central velocity dispersion ctq). The 
measured values of Re, /ie and ao for a sample of E and SO 
galaxies fill a thin plane (with scatter of ~ 0.1 dex) within 



this 3-parameter space instead of spreading out over the 
whole space. The FP can be projected onto any pair of axes 
from the three variables. Examples of these projections are 
the Kormendy relation: fie-^ogRe (Kormendy 1977) and the 
Faber- Jackson relation between luminosity and central ve- 
locity dispersion (Faber & Jackson 1976). 

The distribution of elliptical galaxies in the {^le) — Re 
plane, shown in Figure 1, reveals the existence of two dis- 
tinct families: the "ordinary" family, composed of early-type 
dwarfs, lying in the region defined by 22 {iJ.e{B)} ^ 27 and 
Re ^ 3 kpc and which are always fainter than M{B) = 

— 19.3, and the "bright" family, containing the most lu- 
minous elliptical galaxies (the brightest cluster galaxies - 
BCGs, which host QSOs and Seyfert nuclei), populating a 
relatively narrow strip in the (/Xe) — Re plane, with M{B) < 

— 19.3 and extending over 2 dex in Re- The presence of 
boxy isophotes in these luminous galaxies and the high fre- 
quency of multiple nuclei among BCGs lends support to the 
merging scenario for the "bright" family. The low-luminosity 
disky ellipticals occupy the region between dwarfs and bright 
elliticals in the (fie) — Re plane (—0.5 < logRe < 0.7, 
18 < {fJ-e{B)) < 23). We did not show them in the figure 
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Figure 1. The (fie) — Re plane for different types of objects. Filled triangles and crosses are data from Capaccioli, Caon & D'Onofrio 
(1992). Filled squares are globular clusters from Burstein et al. (1997). Starred symbols are ultracompact dwarf galaxies from Drinkwater 
et al. (2003). The solid line represents the Korincndy relation. Both dashed lines represent the locii of families exhibited in this diagram. 
They schematically represent their overall extension. 



to avoid confusion. Capaccioli, Caon & D'Onofrio (1992) 
interpreted the (/Xe) — Re plane for ellipticals as a logical 
equivalent to the HR diagram for stars, and that galaxies 
might reach the "bright" end of this diagram through suc- 
cessive mergers of low-mass systems. 

The segregation of galaxies in the (ne) — Re plane, to- 
gether with the different behavior of dEs and luminous ellip- 
tical galaxies (Es) in the M{B) — {/le} diagram, where M{B) 
is the absolute magnitude, (see Graham & Guzman 2003), 
was generally interpreted as evidence for distinct formation 
and/or evolutionary processes for the two families of objects. 
However, there is evidence for a continuity, rather than di- 
chotomy, between Es and dEs. According to Graham & Guz- 
man (2003), dE galaxies form a continuous sequence with 
the brighter E galaxies, such that /io (the central surface 
brightness) brightens linearly with M{B) until core forma- 
tion causes the most luminous E galaxies (M{B) < —20.5) 
to deviate from this relation. The different behavior of dE 
and E galaxies in the M{B) — (fie) (or M{B) — He) diagram, 
and the {iie) — Re plane, is unrelated to core formation, and 
is expected from the continuous and linear relation between 
M(B) and no, and M{B) and log(n) (the Sersic index). 

In hierarchical galaxy formation scenarios (as in the 
CDM model), small galaxies (or dark matter halos) are 



the building blocks of more massive galaxies, and should 
have been more numerous in the early Universe. Thus, the 
present-day dwarfs can be seen as survivors of an initially 
much richer population. However, we observe many fewer 
dwarfs today than the predicted number of surviving dark 
matter halos (Klypin ct al. 1999, Moore et al. 1999). There 
also appear to be inconsistencies in the timescales necessary 
to construct large galaxies (Prantzos & Silk 1998), and dif- 
ferences in the stellar populations of large and small galaxies 
(Tolstoy ct al. 2003). These problems point to a need to bet- 
ter understand the physics of dwarf galaxies. 

In the present paper, we concentrate on a single as- 
pect of the complex properties of this family of galaxies, the 
dwarf ellipticals found in the lower left corner of Figure 1 
{{He{B)) ^ 22 and Re kpc). Dwarf ellipticals (dEs) are 
spherical or elliptical in appearance, compact, with high cen- 
tral stellar densities. They are fainter than M{B) = — 18 and 
have low masses {Mtot ^ 10^ ^^0). They are found preferen- 
tially in the vicinity of massive galaxies, usually have little or 
no detectable gas, and are often not rotationally supported. 
When gas is detected, it exhibits an asymmetric distribu- 
tion, is less extended than the underlying stellar component, 
and seems to be kinematically distinct (see Grebel 2001 
and references therein). dEs sometimes show pronounced 
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Table 1. The initial simulation parameters and the relative variations of the total energy and the total angular momentum 
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At N M AE/E% AJ/J% 



Low-mass systems: 



One-component 


ei,2 


= 0.07 


0.025 


Nl,2 


= 90 000 


Ml, 2 


= 9 


0.4 


2.8 


equal mass merging 




















One- component 


ei 


= 0.09 


0.1 


Ni 


= 90000 


Ml 


= 9 


0.7 


2.3 


different mass merging 


€2 


= 0.13 




N2 


= 30 000 


M2 


= 3 






UCDs 


£1,2 


= 0.002 


0.001 


Nl,2 


= 30 000 


Ml,2 


= 0.6 


5.0 


2.0 


Two-component 




= 0.15 


0.15 


Ns 


= 20 000 


Ms 


= 9 


2.5 


1.5 


equal mass merging 




= 0.40 




Nh 


= 60 000 


Mh 


= 27 






Merging of ellipticals 


ei,2 


= 0.07 


0.005 


Nl,2 


= 90 000 


Ml, 2 


= 2 229 


3.0 


1.8 


lying on the KR 




















Merging of objects 


ei,2 


= 0.07 


0.001 


Nl,2 


= 180000 


Mi,2 


= 35 218 


1.3 


4.0 



lying above the KR 



nuclei, with the fraction of nucleated dEs increasing with 
luminosity. The surface density profiles of dEs are best de- 
scribed by Sersic's generalization of de Vaucouleurs R^^'^ 
law and exponential profiles (see Grebel 2001 and references 
therein). The origin of these early- type dwarfs is still un- 
known. Three general clEisses of models have been suggested 
(Grebel 2003): tidal interactions that transform field disk 
galaxies into spheroidal systems; processes associated with 
the birth of small systems (in the CDM model); and frag- 
ments torn from collisions between larger galaxies. However, 
each of these models encounter certain difficulties in repro- 
ducing the observed properties of dE's. 

The effects of dissipationlcss merging on the FP have 
been explored in Capclato, de Carvalho & Carlborg (1995), 
Dantas et al. (2003), and Nipoti, Londrillo & Ciotti (2003), 
among others. In these contributions, the progenitors lay on 
the FP, and the edge-on projection of the FP was found to be 
reproduced by dissipationlcss merging. In contrast to these 
works, our initial models are low-mass early-type galaxies 
which lie neither on the FP nor along the KR. We consider 
the merging of one-component galaxies of both equal and 
different masses, as well as the merging of two-component 
galaxy models, to explore the effects of dark halos on the 
final galaxy characteristics. We also examine mergers of ul- 
tracompact dwarf galaxies (UCDs) , as these constitute a re- 
cently discovered class of very compact objects (Drinkwater 
et al. 2001, Bckki ct al. 2001), populating the FP in a pre- 
viously empty region between the most luminous globular 
clusters and nucleated dwarf galaxies. 

To form a complete picture and compare our results 
with those obtained by other authors, we run a series of 
simulations where the progenitor models lie on the KR. We 
also merge objects which lie above the KR in the {fie) — Re 
plane. These objects are not found in the nearby universe. 



but may have existed at earlier epochs, and are probable 
progenitors for present-day ellipticals. 



2 NUMERICAL METHOD AND INPUT 

2.1 The numerical code 

Our computations utilize an N-body code with a hierarchi- 
cal tree algorithm and multipolc expansion to compute the 
forces, as proposed by Barnes & Hut (1986). The force com- 
putation includes the quadrupole correction terms, following 
Dubinski (1988). 

The set of input parameters (e, the potential softening 
length; A'^, the number of particles; and At, the integration 
timestep) used in the various simulations are listed in Table 
1. Also given in this Table arc the upper limits of the relative 
variation of the total energy, E, and the relative variation of 
the total angular momentum, J, for simulations using each 
set of parameters. The tolerance parameter was the same for 
all simulations, d = 0.8. As discussed in Dantas et al. (2002), 
the choice of the softening parameter e is a compromise be- 
tween spatial resolution and the coUisionless condition of the 
system on evolutionary time scales. It thus depends on the 
particle number A'^, and on the specific density distribution 
profile (see also e.g. Merritt 1996). 

2.2 Galaxy models 

We measure mass and length in units of 10* Mq and 1 kpc, 
respectively. These values, together with the gravitational 
constant G = 1, fix our time and velocity units as 47.2 Myr 
and 20.7 kms~^, respectively. 

The initial models for merging of low- mass galaxies were 
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Table 2. Merging of one-component equal mass models 



Run 


E 


L 


Tper 


20rcr 




logi?e (kpc) 


(Me) 


progemtor(auii j 












0.19-I-/-0.00 


or* r\r\ 1 / r\ r\-\ 

25.00+/-0.U1 


EOl 


-3 





4.2 


144 


275 


0.43+/-0.02 


25.48+/-0.11 


E02 


-7 





4.3 


100 


100 


0.34+/-0.02 


25.00+/-0.10 


EOS 


-5 





5.3 


119 


121 


0.38-I-/-0.02 


25.23+/-0.12 


E04 


-3 


1 


24.3 


140 


150 


0.43+/-0.02 


25.50+/-0.11 


E05 


-5 


1 


11.3 


117 


120 


0.38-I-/-0.03 


25.21+/-0.13 


E06 


-7 


1 


6.8 


102 


120 


0.34-I-/-0.02 


25.02-I-/-0.10 


E07 


-3 


2 


24.3 


140 


150 


0.43+/-0.02 


25.46+/-0.11 


EOS 


-5 


2 


11.6 


91 


150 


0.38+/-0.02 


25.20+/-0.10 


E09 


-7 


2 


6.8 


102 


175 


0.34-I-/-0.02 


25.00+/-0.10 


ElO 


-1 


3 


126.3 


174 


200 


0.47+/-0.02 


25.70+/-0.12 


Ell 





1 


2.6 


166 


170 


0.49+/-U.03 


25.87+/-0.15 


E12 


0.5 


1 


2.5 


178 


182 


0.49-I-/-0.03 


25.87+/-0.15 


IljLo 


i 


1 
i 


K 

Z.O 


101 


zUU 


n 1 / n hq 
U.OU-h/-U.Uo 


OK OK 1 / n 1 /I 

zo.yo-h/-u.i4 


EOI2 


-3 





5.8 


193 


234 


0.62-I-/-0.03 


26.76+7-0.16 


E022 


-7 





4.2 


82 


150 


0.47+/-0.05 


24.93+/-0.24 


E052 


-5 


1 


11.6 


109 


150 


0.54+/-0.04 


25.27+/-0.20 


EO62 


-7 


1 


7.0 


83 


125 


0.46+/-0.04 


24.87+/-0.18 


EO82 


-5 


2 


11.6 


96 


231 


0.55+/-0.04 


25.33+/-0.20 


E092 


-7 


2 


7.0 


83 


150 


0.47+/-0.04 


24.90+/-0.19 


EOI3 


-3 





8.9 


416 


546 


0.83-I-/-0.03 


26.19-I-/-0.17 


EO23 


-7 





4.2 


85 


150 


0.59+/-0.05 


24.78-I-/-0.26 


E053 


-5 


1 


11.6 


84 


200 


0.71-I-/-0.04 


25.37+/-0.22 


E06:5 


-7 


1 


7.0 


134 


359 


0.58+/-0.05 


24.71+/-0.24 


E093 


-7 


2 


7.0 


87 


300 


0.59+/-0.06 


24.75+/-0.27 



Table 3. Merging of one-component different mass progenitors 



Run 


E 


L 




20Tcr 




logiie (kpc) 


(Me) 


progcnitorl(a!tyi ) 












0.19+/-0.00 


25.00+/-0.01 


progemtor2 (dwi ) 












0.18+/-0.00 


26.16+/-0.01 


DOl 


-3 


1 


35.0 


110 


164 


0.36+/-0.02 


25.60+/-0.10 


D02 


-5 


1 


16.3 


99 


152 


0.34-I-/-0.02 


25.46-I-/-0.09 


D03 


-7 


1 


9.9 


91 


137 


0.30-I-/-0.01 


25.28+/-0.07 


D04 


-3 


2 


35.0 


110 


164 


0.36-I-/-0.02 


25.58+/-0.10 


D05 


-5 


2 


16.3 


100 


152 


0.33+/-0.02 


25.44+/-0.09 


D06 


-7 


2 


9.8 


92 


140 


0.30+/-0.02 


25.28+/-0.08 



chosen from the region in Figure 1 occupied by ordinary el- 
lipticals and early- type dwarfs (the "ordinary" group) . Both 
one- and two-component models were constructed as de- 
scribed below. 



The initial one-component galaxies were modelled us- 
ing the potential-density pair for spherical galaxies given in 
Hernquist (1990) : 



$(r) = 



GM 
r -|- a 
M 



(1) 

(2) 



r (r + a)'' ' 

where M is the galaxy mass and a is the scale length. 

In the case of two-component models (the luminous 
galaxy and its dark matter halo), the halo density profile 
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Table 4. Merging of ultracompact dwarf galaxies 



Run E L Tper 20Tcr Tend logiJe (kpc) (^te) 



progenitor(Mcd) -1.65+/-0.00 18.75+/-0.01 



FOl 


-1 


1 


0.84 


0.96 


2.0 


-1.32+/-0.03 


19.69+/-0.13 


FOI2 


-1 


1 


1.80 


1.80 


4.0 


-1.03+/-0.04 


20.44+/-0.21 


FOI3 


-1 


1 


3.22 


3.55 


7.3 


-0.74+/-0.05 


21.14+/-0.23 


FOI4 


-1 


1 


5.78 


6.82 


18.5 


-0.47+/-0.04 


21.80+/-0.22 



Table 5. Merging of two equal mass progenitors with dark halo 



Run 


E 


L 


Tper 


20Tcr 




logi?e(kpc) 






progenitor 












0.19+/-0.00 


25.01- 


h/-0.01 


HOI 


-3 


1 


55.8 


245 


700 


0.37+/-0.03 


25.34- 


h/-0.17 


H02 


-5 


1 


26.0 


204 


700 


0.36+/-0.04 


25.31- 


h/-0.17 


H03 


-7 


1 


15.7 


174 


600 


0.35+/-0.04 


25.26- 


h/-0.18 



was represented as a truncated isothermal sphere (Hernquist 
1993): 

Here Mh is the halo mass, 7 is the core radius, and Vc is the 
cutoff radius. The normalization constant a is defined by 

a = (1 - V¥gea;p(g^)[l - er/(g)])"\ (4) 

where q — 7/rc. The cumulative mass profile, Mh{r), and 
potential, ^(r), corresponding to p{r) are 



in dwarf galaxies, these models are consistent with current 
ideas on the structure of these objects. 

The effective radii of UCDs range from 0.015 kpc 
(logiie = -1.82) to 0.03 kpc {logRe = -1.52) and their 
mean surface brightness within Re ranges from 18 to 19. 
The UCD radial profiles are well described by de Vau- 
couleurs profiles (Drinkwater et al. 2003). For the initial 
UCD models, we used a Hernquist sphere with Mucd ~ 0.6 
and aucd = 0.014. 

Finally, for the progenitor models lying on and above 
the KR, we constructed models with Mkr = 2229, and 
Makr ~ 35218, respectively, with both models having a = 
1. 



Mn{r) = ^^ / Jf^^dx, 5 

^ J 



. , GMh{r) GMh{r)a / , n2 21 /«n 

$h(r) = ^ -I Ei[-{r re) -g , (6) 

r y/nrc 

where Ei{z) is an exponential integral. 

We constructed single-component models representing 
early- type dwarfs galaxies with total masses of Mdwi ~ 9 
and Mdw2 = 3, both having the same scale length a — 1 
{Re X 1.8a). 

The two-component model for early-type dwarfs con- 
sists of a luminous component identical to that for one- 
component dwi model described above, with the dark mat- 
ter halo more massive {Mh/Ms — 3) and less concentrated 
(7/a = 3) than the stellar component. The dark compo- 
nent parameters for these models are Mh ~ 27, 7 = 3, and 
Vc — 12. Given the lack of observational data on dark matter 



2.3 Initial conditions 

The encounter of two non-rotating spherical galaxies may 
be characterized by the dimensionless energy and angular 
momentum: 



L = , (8) 

where (v^) = \/{vl){vl), = y^nun^. Here (u^) is the 
internal mean-square velocity, rn is t he half-mass radius of 
a galaxy (Binney & Tremaine I^S^))- The indices denote 
each of the initial galaxies and p, is the reduced mass of the 
system. 

The time required for two galaxies to merge is a function 
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Table 6. Merging of ellipticals lying on the KR 



Run E 


L 


20Tcr 




logiJe (kpc) 


(Me) 


progenitor(A'_R) 








0.18+/-0.00 


19.00+/-0.01 


POl -3 


1 


7 


18 


0.44+/-0.03 


19.54+/-0.13 


P02 -5 


1 


6 


18 


0.39+/-0.02 


19.30+/-0.11 


P03 -7 


1 


5 


18 


0.35+/-0.02 


19.08+/-0.11 


POI2 -3 


1 


11 


22 


0.66+/-0.04 


19.92+/-0.20 


PO22 -5 


1 


8 


24 


0.56+/-0.04 


19.41+/-0.19 


P032 -7 


1 


6 


20 


0.49+/-0.04 


19.01+/-0.21 


. Merging of ellipticals lying above the KR in 


the (He) - Re 


plane 






Run E 


L 


20 Tcr- 




logRe(kpc) 


<Me> 



progenitor(AiirR) 



A04 
A05 
A06 
A07 
AOS 
A09 



-3 
-5 
-7 
-3 
-5 
-7 



1 
1 
1 
2 
2 
2 



1.74 
1.46 
1.42 
1.74 
1.46 
1.26 



0.19+/-0.00 16.03+/-0.00 



0.44+/-0.02 
0.39+/-0.02 
0.34+/-0.02 
0.44+/-0.02 
0.38+/-0.02 
0.34+/-0.02 



16.55+/-0.12 
16.28+/-0.11 
16.05+/-0.10 
16.53+/-0.11 
16.26+/-0.11 
16.04+/-0.11 



of the initial position of the binary orbit in the {E,L) plane 
defined by equations (7) and (8) (Binney & Tremaine 1987). 
We considered only those pairs of {E, L) which correspond 
to rapid mergers (less than a Hubble time). 

To completely determine the orbital elements of the en- 
counter, we must define a third parameter: 



A. 



2GM 



(9) 



which depends only on the dynamical structure of the ini- 
tial galaxies. The initial separations and velocities for the 
encounters are selected from a grid of E, L, A values. The 
initial separation of the models was taken at the apocenter 
position for bound orbits and ~ 4r/t for unbound ones. 

We follow the dynamical evolution of each merger until 
the resulting system is virializcd, which typically occurs on 
a time scale shorter than 20 Tcr after the first encounter 
between two galaxies, where Tcr is the crossing time of the 
final object (Tcr- = GM^^^ /{2Ef/^). The values for this 
timescale are given in Tables 2-7, column (5). 

Tables 2-7 summarize the initial conditions and charac- 
teristic parameters of the simulations. In Table 2, the merg- 
ers of one-component equal mass galaxies are labelled as 
follows: 

• EOl - E13: first- generation mergers: simulations of the 



encounter of two single- component equal mass progenitors 
with different pairs of {E, L); 

• EOI2 - E092: second-generation mergers: simulations of 
the encounter of two identical end-products of first genera- 
tion mergers from the previous set of simulations. Progen- 
itors are selected such that E0i2 = (EOi -|- EOi), with i = 
{1,2,5,6,8,9} ; 

• EOI3 - EO93: third- generation mergers: encounters be- 
tween two identical second generation mergers with the pro- 
genitors selected such that EOis = E0i2 -|- E0i2, with i = 
{1,2,5,6,9}. 

The other simulations (Tables 3-7) with similar sub- 
scripts but a different starting letter, are as follows: 

• D - mergers of different mass one-component mass 

galaxies; 

• F - mergers of UCD models; 

• H - mergers of two-component equal mass galEixies; 

• P - mergers of one-component models lying on the KR; 

• A - mergers of one component models lying above the 
KR. 

In Table 6, POI2 - P032 are the second generation merg- 
ers which, as for the E mergers, had progenitors selected 
such that P0i2 = (POi + POi), with i = {1,2,3}. Similarly, 
in Table 4, FOlj = (FOlj-i -\- FOlj-i), where j = {2,3,4}. 

Columns (2) and (3) in Tables 2-7 give the initial di- 
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Figure 2. Merging of equal mass progenitors. * - the initial pro- 
genitor. Circles - the 1st generation progenitors; squares - the 
2nd generation progenitors; triangles - the 3rd generation pro- 
genitors. The dashed line corresponds to the line in Fig.l which 
encompasses the region of bright ellipticals below and it is par- 
allel to the KR. The rectangular area embraces the locus of the 
dwarf galaxies (see Fig.l). 
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Figure 4. Same as in Figure 2 but for merging of ultracompact 
dwarf galaxies. 
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Figure 5. Same as in Figure 2 but for merging of two equal mass 
progenitors with dark halo. 



Figure 3. Same as in Figure 2 but for merging of two different 
mass progenitors. 



mensionless orbital energy and corresponding angular mo- 
mentum. Columns (4) and (6) give Tper, the predicted time 
interval before the two-body pericenter, and Tend, the total 
time elapsed to the end of the simulation, respectively. 



3 RESULTS 

We followed the procedure described in Capclato, dc Car- 
valho & Carlberg (1995) to quantify the observables involved 
in the KR: the effective radius Re, containing half of the total 
projected mass of the system, and the mean surface density 



within Re, Ee = M(< Re)/'KRl, so that (fie) = -2.51ogEe. 
These quantities were estimated as the median of 500 ran- 
dom projections. The errors quoted for each parameter are 
the rms based on the quartiles of the distribution defined by 
the 500 random projections. Columns (7) and (8) in Tables 
2-7 list Re the effective radius, and (yUe), the mean surface 
brightness inside Re, for the merger remnant of each simu- 
lation. These quantities are shown in Figures 2-7. 

In Table 2, the values of E and L defining the encounter 
span the region of rapid merging in the E — L diagram (Bin- 
ney & Tremaine 1987). Although arbitrary, this choice allows 
us to examine a particular scenario where mergers occur over 
a short timescale. Some global trends can be seen from this 
set of simulations. First, (//e) always increases (lower sur- 
face brightness) with Re, and the rate of increase depends 
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Figure 6. Merging of ellipticals lying on the KR. * - the initial 

progenitor. Circles - the 1st generation progenitors; squares - the 
second generation progenitors. The dashed lines correspond to 
those in Fig.l. 
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Figure 7. Merging of ellipticals lying above the KR in the {fie) — 
Re plane. * - the initial progenitor. Circles - the 1st generation 
progenitors. The solid line represents the KR. The dashed lines 
correspond to those in Fig.l. 



on the initial orbital energy. Second, the presence of angular 
momentum in the initial conditions does not significantly 
change the properties of the merger remnants, a trend al- 
ready noted by Capolato, do Carvalho, & Carlborg (1995). 
Therefore, we restrict ourselves to a smaller range of L for 
each E value in the further simulations described here. We 
use the same {E, L) for the simulations presented in Tables 
3-7 as for those in Table 2, except for the E = -1 case in- 
cluded in Tabic 2, for which Tper is very large compared to 
the other simulations. 

In the simulations presented in Tables 3-7 we restrict 
ourselves only to the first generation (except for the merging 



of objects lying on the KR - Table 6, and the merging of 
UCD models - Table 4), since the objects produced in the 
first generation tend always to move in the same direction 
in the (^e) — Re plane as the first generation systems listed 
in Tabic 2. Judging by the behavior of the the second and 
third generation families displayed in Figure 2, we can safely 
assume that further merging will also move the end-products 
to the bottom right in the {jie) — Re plane. 



4 DISCUSSION 

In Capolato, dc Carvalho 



; Carlborg (1995), 

-ill : 



the progen- 



have used low- mass progenitors ranging from 0.6x10® to 
9xlO®M0, since the main goal was to study whether or 
not the low-mass systems also move along the FP. How- 
ever, the loci of the progenitors studied hero arc not along 
the KR (a projection of the FP). From our simulations, we 
see that present-day ellipticals cannot be formed by merging 
present-day dwarfs, unless a considerable amount of dissipa- 
tion is involved. This result, together with other inconsisten- 
cies present in the properties of dwarfs and bright galaxies 
(see Tosi 2003), makes it nearly impossible that the major 
building blocks of present-day ellipticals had properties sim- 
ilar to those of dwarfs galaxies observed today. 

The global trend observed in Figures 2-7 is that the 
end-products move down and right in the (/Xe) — Re plane. 
On one hand, this shows that the only way to merge dwarfs 
and obtain bright ellipticals is by invoking a large amount 
of dissipation. On the other hand, the same trend suggests 
the possibility of making dwarfs by merging UCDs (clearly 
shown in Figure 4 and Figure 8), as well as the possible 
formation of UCDs and dwarfs by merging globular clusters. 
Fellhauer & Kroupa (2002) have also shown that objects 
such as UCDs can be formed by merging globular clusters. 
Interestingly, in the dissipationless simulations by Dantas et 
al. (2003), the line going through the positions of the mergers 
of all of their three generations with —4 < E < -1-0.5 has 
the same slope as the KR. However, their simulations do 
not reconstruct the KR entirely. In essence, all these results 
show how dissipationless merging moves the systems in the 
(/ie) — Re plane, regardless of their mass. 

The structural parameters of the end-products seem to 
strongly depend on the initial orbital parameters {E, L, and 
A). A recent study done by Khochfar & Burkert (2003) on 
the orbital parameters of major mergers of CDM halos in- 
dicates that most of the encounters are nearly parabolic. 
However, as shown by simulations Ell, E12, and E13, which 
correspond to parabolic and hyperbolic orbits, the trend is 
the same as observed for the other cases, namely, the end- 
products still move down and right in the (ne) — Re plane. 
The global trend is independent of the E range, no matter 
what kind of orbits we consider. 

The simulations where dark halos were added to the 
low-mass spherical galaxies show that the merger remnants 
also shift down-right in the (/Xe) — Re plane (Table 5 and 
Figure 5) , indicating tliat the presence of dark halos does not 
make any appreciable difference as far as the global trend 
we find is concerned. 

We also simulated the merging of a pair of objects ly- 
ing on the KR. The end-products also lay on the KR within 
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Figure 8. The evolution of dwarf galaxies in the {fie) — Re plane. Merging of ultracompact dwarf galaxies is shown by large arrows. 
Different symbols represent the evolution of a dwarf galaxy due to the effects of dissipationless encounters for three subsequent merger 
generations. The solid line represents the Kormendy relation. Both long-dashed lines correspond to those in Fig. 1 which define the lower 
and the upper boundaries for bright ellipticals. 



the scatter of the observed KR, and taking into account the 
simulation errors (Table 6 and Figure 6), are in agreement 
with the results obtained by Capelato, do Carvalho & Carl- 
berg (1995), and Dantas et al. (2003). The simulations where 
the progenitors were above the KR relation yielded the same 
overall result, moving down and to the right in the {fj,e} — Re 
plane. 

In Figure 8 we show with different symbols the evo- 
lution of a progenitor duo to the effects of dissipationless 
encounters (see Table 2 and Figure 2) for three subsequent 
merger generations. The resultant vectors point to the locus 
of the bright ellipticals defining the KR. From the scheme 
presented in Figure 8 we see that we cannot obtain a bright 
elliptical by merging low surface brightness ((/ie) = 25 — 26) 
dwarfs. However, it is possible to reach the KR for bright 
ellipticals by merging higher surface brightness dwarf galax- 
ies ((pe) = 22 — 23) in a few (two or three) merger episodes, 
in agreement with Bower, Kodama & Terlevich (1998). We 
can imagine that at the initial merging stages, the dwarfs' 
main interaction process is due to dissipational effects. In 
the (/^e) — Re diagram, this process is primarily represented 
by the dwarfs' upward movement. Then, when these dwarfs 



reach the upper part of the dwarf locus, they will move ac- 
cording to the scheme in Figure 8. 

There arc some caveats that should be noted. Dissipa- 
tion could have a significant impact on the final position of 
an end-product. Also, it is important to bear in mind that 
our simulations start from equilibrium systems representing 
a realization of a system which was previously formed. More- 
over, the present merger rate differs from the one which was 
in the past. 

We plan on doing more extensive simulations using dif- 
ferent models for the initial low-mass ellipticals (e.g. the ex- 
ponential model or the more general Sersic's law) as well as 
those for the dark halo. To study how the initial spin of our 
models could affect the end-products is also an important 
part of the future development. 
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